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Abstract 


A feedback control system has been designed and simulated for transverse 
injection into a laminar, supersonic flow over a flat plate. This control system 
employs the injection total pressure of nozzle as a control input to the system. Two- 
dimensional Navier-Stokes equations are solved with an explicit scheme to simulate 
the flow field. Numerical results are presented for different desired static pressure 
ratios and with various values of controller constants. Robustness of the controller 
with respect to the ffeestream Mach number is demonstrated. A transfer function of 
the linearized system is obtained from the numerical results. 
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Chapter 1 


Introduction 

1.1 Literature Survey 

The significance of studying the flow field set up by normal injection of secondary 
gas into a supersonic primary flow has been a problem of interest because of its 
application in fuel injection into a supersonic combustor, thrust vector control of 
rocket motors and in high speed flight vehicle reaction control jets. A number of 
experimental studies were carried out initially to understand the complex flow field 
set up by normal injection into supersonic flow. Analytical and numerical studies 
were also been carried out to describe this complex flow field. This section gives a 
brief note on the various studies conducted and the advancements achieved in 
understanding this flow field. 

Dorrance* studied the use of introducing a secondary fluid into a primary flow, 
which has a hot gas boundary layer for cooling the surface. He developed equations to 
study the effect of mass transfer and chemical reactions upon the laminar boundary 
layer properties by solving the boundary layer equations along with some chemical 
equations. But his theory is basically to study the effect of introducing the secondary 
fluid for cooling purposes. 
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The first experimental investigation to study the flow field set up by normal 
injection into supersonic flow was carried out by Spaid and Zukoski They 
conducted a series of experiments at fireestream Mach numbers of 1.38 to 4.54 with 
gaseous nitrogen, argon and helium as injectants and proposed a scaling parameter 
which allows for the simple prediction of side forces generated by transverse 
injection. Schetz and Billing"^ studied the behaviour of gaseous jets issuing 
transversely into a supersonic primary stream analytically and the analysis was 
applied to a configuration study of a combustion chamber for a hypersonic ramjet. 

George Lee ^ used the boundary layer theory along with hypersonic small 
disturbance theory to derive the formulas for mass-transfer distribution, pressures, 
heat transfer, skin friction and boundary layer displacement thickness for the 
hypersonic flow over a blunt flat plate with surface mass transfer. He proved that the 
surface pressures would be increased by mass transfer and with higher ratio of wall to 
stream temperature. Mass transfer will strongly influence the heat transfer and makes 
it to approach zero for high rates of mass transfer. Spaid and Zukoski ^ conducted flat- 
plate experiments with normal sonic jets at external flow Mach numbers of 2.61, 3.5 
and 4.54 and developed an analytical model of the flow field by applying momentum 
conservation to a control volume at the jet nozzle exit. They correlated the data 
obtained from experiments with finite-span slots and demonstrated that the effective 
jet penetration height and the slot span are the important characteristic dimensions of 
such flowfields. Spaid ’ conducted experiments with normal, sonic jets and forward- 
facing steps at external flow Mach numbers of 2.5 to 13, and Reynolds number based 
on running length of 7.5 x 10^ to 5.5 x 10®. 

The flow field was simulated numerically using different turbulence 
models and various forms of Navier-Stokes equations. Napolitano ® used the triple- 
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deck equations to solve the problem of strong slot injection into a supersonic laminar 
boundary layer and solved them for a wide range of wall injection velocities. Weidner 
and Drummond ^ solved the full two-dimensional Navier-Stokes and species 
equations to analyze the turbulent reacting flow field in a two-dimensional scramjet 
engine configurator, especially in engine inlet and combustor. They also conducted 
a parametric study of flowfields near staged perpendicular injectors by solving the 
two-dimensional elliptic Navier-Stokes equations. They used MacCormack time-split 
algorithm with algebraic eddy-viscosity model along with a chemistry model to study 
the mixing phenomena of hydrogen fuel and air. They suggested the use of staged fuel 
injection for scramjet engine. 

Rizetta ' ‘ numerically simulated the slot injection experiments of Aso et al 
at freestream Mach number of 3.7 with turbulence modeling simulated by eddy 
viscosity model. He assessed the adequacy of the turbulence model with regard to 
nozzle width and slot pressure ratio for jet injection flowfields. Chenault and Beran 
studied the feasibility of compressible k-8 and second-order Reynolds-stress 
turbulence model for transverse injection into a supersonic flow by integrating the 
Favre-averaged Navier-Stokes equations. Results of his study indicated that the 
Reynolds-stress turbulence model correctly predicts mean flow conditions and the k-s 
results are less consistent with the experimental results than those associated with the 
Reynolds-stress turbulence model. But the vorticity phenomena upstream of the jet 
were best observed with k-e model. 

While most of the above discussion is related to transverse injection into 
supersonic flow, there are many other applications where mass injection into external 
flow will be done through porous surfaces or by vectored or tangential injection. 
These types of injection finds application in a number of important boundary layer 
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flow problems such as film and transpiration cooling of rocket engines, turbine 
blades, surfaces of high-speed aerodynamic bodies and for skin friction reduction and 
boundary layer control on high-speed aircraft. The following paragraphs will give a 
brief note on the studies conducted in these t>pes of flows. 

Kenworthy and Schetz studied the use of injection of a gas through a wall 
slot parallel to a supersonic stream for skin-friction reduction on high-speed aircraft. 
Inger and Swean conducted a theoretical investigation of self-similar isobaric 
laminar two-dimensional boundary layer flows for a wide range of conditions for 
suction, as well as injection, including the effects of heat transfer. They presented 
detailed results for skin friction, velocity and enthalpy profiles and momentum and 
displacement thickness integral properties. Nilson and Tsuei solved the boundary 
layer equations numerically by a finite difference method to study the film cooling 
effectiveness for oblique injection of coolant through single or multiple wall slots into 
a high-speed laminar flow. Schetz and vanOvereem made an experimental study to 
compare the effectiveness of tangential slot injection and porous wall injection for 
skin friction reduction. Inger and Swean extended their work of vectored 
injection and suction to study the effects of self-induced pressure gradient caused by 
the boundary layer displacement thickness growth. They presented the results 
showing the effects of both injection velocity and vector angle on a variety of physical 
properties. Hefner for the first time, conducted experiments with two different slot 

geometries to study whether any simple modifications to the slot configuration can 

20 

produce improved cooling effectiveness and skin friction reduction. Brandies 
numerically studied the flow field within a finite, two-dimensional slot through which 
fluid can either be added or subtracted and also studied the effects of the size of the 
slot on the flow field. Connor and Sheikh conducted a computational fluid 
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dynamics study which focuses on injecting a heated film of air through a rearward 
facing slot parallel to a flat plate in a Mach 3 divergent nozzle. 

Dershin et al conducted an experimental investigation to calculate the skin 
friction on a porous plate with mass injection at a supersonic Mach number of 3.2 and 
with turbulent Reynolds number. But, the concentration was more on measurement of 
skin friction and predicted the skin friction reduction with mass injection. Thomas 
used an inviscid, vortical flow model for analyzing the compressible flow over a finite 
porous plate in supersonic stream, with massive injection at a uniform rate over the 
surface. The above finds its application in cooling the surface during re-entry from an 
interplanetary mission. Clark et al conducted an experimental investigation to 
determine the local and downstream cooling effectiveness associated with injection of 
water through slotted and porous surfaces into a turbulent, supersonic boundary layer. 
El-Mistikawy and Werle developed a stable second-order accurate finite-difference 
scheme to provide solutions for boundary layer flows in the presence of massive 
injection through a porous surface. 

The methods described above, i.e., transverse flow injection, parallel and 
vectored injection and mass injection from porous surface come under the category of 
active flow control where some additional input is used to get the desired effect which 
can be either skin-friction reduction, thrust vector control, fuel injection into 
supersonic combustion engine, cooling surfaces of high-speed aircraft etc. So, the 
interest in the active flow control for the above applications has stimulated the recent 
development of innovative actuator designs that create localized disturbances in a 
flow field. The following paragraphs will give a brief not on the recent developments 
achieved in actuator designs and modem flow control methods. 
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Nae presented results for numerical simulation of synthetic jet actuator and 
the flow around the NACA 0012 airfoil for a given control law and for circular 
cylinder using active flow control. Calkins and Mabe developed a synthetic jet 
actuator to test in a 1/10 scale V22 2D wind tunnel model and in a 1/10 scale 3D 
powered model. S 3 mthetic jets result from oscillating diaphragm in an enclosed 
space, having small orifices at the top. The jet will be synthesized from the working 
fluid of the flow system requiring no external fluid source, thus they are commonly 
referred as zero-net-mass-transfer-flux actuators. They can be controlled electro 
statically or using piezo electric materials with frequencies in the range of 0.5-20 
ICHz. These have been found to be more effective than other traditional mechanical 
control systems for flow control applications. Another means of flow control can be 
achieved by pulsed injection. Vakili et al have experimentally and numerically 
investigated the use of such a pulsed injection in a nozzle. Their objective was to 
increase the fluidic blockage produced by a transverse injector into the exhaust stream 
near the throat. Their experimental results depicted that the jet penetration was 
increased with pulsation of the jet. Also the results for various exit geometry jets 
indicated an improvement in the effective blockage for periodic injection compared to 
steady state injection. 

The technique of flow control by magnetogasdynamic (MGD) forces was 
simulated by Gaitonde and Poggie In this, a plasma which may have been created 
naturally by the presence of shocks and/or artificially through some form of energy 
input like laser beam into the flow will be used for flow control. They developed a 
simulation tool by integrating the Maxwell and compressible Navier-Stokes equations 
with the assumptions of non-ideal magnetogasdynamics and investigated various 
issues concerning the numerical simulation of such a hypersonic flow control 
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problem. Lmeberry et al studied the prospects of active magneto hydrodynamic 
(MHD) boundary layer control and leading edge heat transfer management with 
MHD. They also qualitatively described MHD flow control inlet configuration 
concept and its operation. Also an active flow control by electro hydrodynamic 
(EHD) actuator was studied by Leger et al for low velocity airflow along a flat 
plate. Such an actuator can be used to accelerate the airflow tangentially and very 
close to the wall, in order to modify the airflow profile within the boundary layer. Its 
advantage is that it directly converts electrical energy into kinetic energy without 
moving any mechanical part. Its response time is very short and enables a control as a 
function of time at high frequency. Its disadvantage is its low efficiency of energy 
conversion and also generation of extraneous gases such as ozone in air. Riggins and 
Nelson numerically studied a flow control concept in which energy will be 
deposited at a specific point by solving the full Navier-Stokes equations. They 
conducted study on centerline and off-centerline power deposition at a point upstream 
of the two-dimensional blunt body at Mach 6.5 at 30 km. altitude. This flow control 
concept can be effectively used to improve the performance and to stabilize and 
control hypersonic vehicles. 

While many experimental and numerical studies of active flow control 
described above have been carried out, an active flow control problem of transverse 
injection with a feedback has not yet been reported. 

1.2 Transverse Injection into Supersonic Flow: Flow field 
Description 

A two-dimensional flow field set up by injection of a secondary fluid normal to an 
external supersonic stream is shown in Fig. (l-l)- In this example, it is assumed that 
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the jet is sonic, underexpanded and the obstruction to the external flow formed by the 
jet is larger than the undisturbed boundary layer thickness. Expansion of the jet occurs 
rapidly as it leaves the nozzle. The interaction of the secondary jet causes a strong 
bow shock to form upstream of the jet which provokes the boundary layer to separate. 
Another shock wave will be generated at the separation point, called as separation 
shock, which intersects the jet induced bow shock. Upon exiting the nozzle, the jet is 
decelerated by a normal shock, often referred as Mach disk, to obtain the pressure 
equilibrium. Downstream of the slot, the boundary layer reattaches resulting in a 
recompression shock and a subsequent separation region. 


separation shock 



boundary 

layer 


upstream vortices 


downstream vortices 


Figurel.l: Schematic of the flow field 


The surface static pressure distribution for transverse injection'^ is shown in 
Fig. (1.2). The static pressure rises because of boundary layer separation (region 1), 
followed by a plateau (region 2) because of constant flow conditions of the upstream 


8 



vortices, and rises again (region 3) in the immediate vicinity of the jet because of the 
bow shock. Immediately downstream of the jet is a region (region 4) where both the 
flows will separate and pressure falls below the ambient value. This is followed by a 
pressure rise (region 5) because of recompression shock and boundary layer 
reattachment. 



Figure 1 .2: Schematic of Surface Pressure Distribution 

Spaid and Zukoski ^ reported that when the boundary layer is laminar, the 
separation region is large and the angle between the separated flow and the wall will 
be very small. For the case of a turbulent boundary layer, the separation region will be 
small and the shock produced by separation is usually sufficiently strong to be 
observed. 
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1.3 Passive and Active Control 


Passive control is a control strategy employed for a system where external inputs are 
not required for the control of the system. The system may be consisting of 
subsystems, which effectively perform the control function, without the need for 
distinct components for sensing and actuation. Passive flow control is normally 
achieved by incorporating changes in the flow geometry (e.g. introduction of fixed 
slots, flaps or sharp edges in the flow field) 

Apart from passive control, there are control strategies which require distinct 
components, including actuators and sensors. This strategy is called as ‘active 
control’. This strategy uses additional inputs for the control of the system. An active 
control can be either an open loop or closed loop. An open loop active control system 
is shown in Fig. (1.3). 



CONTROLLED 



SYSTEM 


f 

(PLANT) 

y 


Figure 1.3: Active Control (Open loop) 


The output of the above system may not be a desired output. If the controlled system’s 
behavior is precisely known, then it is possible to achieve the desired performance 
through open-loop control. However, this requirement is rarely met. 

Closed loop (or feed back) control is desirable if the controlled system’s 
characteristics are unsatisfactory (such as instability, low-damping, large overshoot. 
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steady-state error etc.), and the closed-loop system (if properly designed) can improve 
the overall performance. Usually, a closed-loop system employs sensors, a feedback 
path, a forward path and actuators as shown in Fig. (1.4) 


Forward path 



Feedback path 


Figure 1 .4: Closed loop Active control 

Here, f(t) is the desired output vector, y(t) the actual output vector, z(t) the sensed 
output vector and u(t) the control input vector. The design of the controller must be 
carefully carried out to achieve the desired closed-loop performance. Modem control 
theory^^'^^ is based upon either a transfer function or a state-space description of the 
controlled system (plant), and is largely restricted to linear systems. However, most 
non-linear systems can be linearized about their equilibrium points, and thus the 
concepts of modem control theory can be applied to design a controller for non-linear 
systems where only small deviations from an equilibrium point are expected. This is 
the approach adopted in the present thesis, where the plant is described by nonlinear 
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partial differential equations (Navier-Stokes equations) and linearized transfer 
function approximation is used for designing the controller. 


1.4 Types of Feedback Controllers 

The commonly used feedback controllers for SISO plants are 

1 . Proportional controller 

2. Integral controller 

3. Proportional-plus-integral controller 

4. Derivative controller 

5. Proportional-plus-Derivative controller 

6 Proportional-plus-Integral-plus-Derivative controller 


1.4.1 PID controller 

The transfer function of a PID controller is of the form 


Ufe) 

E(s) 


Kp+ K,s + 


K; 


( 1 . 1 ) 


and its block diagram is shown in Fig. (1.5) 



U(s) 

. . . 'w-. 

G(s) _ 

^ 



Y(s) 


Figure 1.5: Block Diagram of a PID Controller 
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where U(s) = X { u(t) }, E(s) = X { e(t) }, Y(s) = X { e(t) } 

u(t) is the control input, y(t) is the output and e(t) is the error signal. G(s) is the plant’s 
transfer function and Kp, IQ and K, are constants. The integral term, K,/s, introduces a 
pole at s=0, thus increasing the type of the closed loop system by one. This reduces 
the steady state error of the original system. The proportional term Kp, and derivative 
term, BQ s, can be used to improve the transient response of the system by placing two 
zeros at suitable locations. They can be made to yield a system with fast rise time and 
with little or no overshoot because of increased system damping. Thus a PID 
controller will have the advantages of all the three controllers and thus can be used to 
improve both the transient and steady state response of the system. 
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Chapter 2 


Active Control Strategy and Governing Equations 

2.1 Formulation of Navier-Stokes Equations 

2.1.1 Continuity Equation 

The conservation of mass law applied to a perfect gas passing through an 
infinitesimal, fixed control volume yields the following equation of continuity 

^ + V-(/?V) = 0 (2.1) 

ot 

The first term in this equation represents the rate of increase of the density in the 
control volume and the second term represents the rate of mass flux passing out of the 
control surface (which surrounds the control volume) per unit volume. 

For a Cartesian coordinate system, where u, v, w represent the x, y, z 
components of the velocity vector V, the above equation becomes 



2 . 1.2 Momentum Equation 

Newton’s second law applied to a perfect gas passing through an infinitesimal, fixed 
control volume yields the following momentum equation 


DV . ^ 9 

p = pf-Vp + — 

Dt ^ dx. 


M 


du, dUj 
V y 


2 du^ 

3 ^ dx. 


or 




(2.3) 


(2.4) 


where is the viscous stress tensor and 5 is the Kxonecker delta function {5 =1 


if i = j and <5',^ = 0 if i j ); ui, U 2 , U 3 represent the three components of the velocity 

vector V; xi, X 2 , X 3 represent the three components of the position vector. 

The first term on the right-hand side of the above equation is the body force 
per unit volume. Body forces act at a distance and apply to the entire mass of the 
fluid. The most common body force is the gravitational force. In this case, the force 
per unit mass ( f ) equals the acceleration due to gravity vector g 

pi = pg (2.5) 

The second and third term on the right-hand side represents the surface forces per unit 
volume. These forces are applied by the external stresses on the fluid element. The 
stresses consist of normal stresses and shearing stresses. 

For a Cartesian coordinate system, the Navier-Stokes equations can be written 
in conservation-law form as 

+ + ^(puw-Tj=pf^ 

dtdx oy oz 

dt OX oy oz 

8(^+A(^,„-r„)+A(pvw-r„)+ +f.-r„)=p/, 

dtdx ' oy oz 
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where the components of the viscous stress tensor are given by 






T = — U. 
zz 


du 

dv 


dx 

dy 

dz ^ 

^ dv 

du 

dw^ 


dx 

dz 

> 

dw 

du 

dv 

'^dz 

dx 

dy 


P- 


du dv 
dy dx 


= T.. 


dw du 
dx dz 


= r. 


y^ 




dv dw 
dz dy 


= r. 


zy 


and 


f = /.i + 4 j + /.k 

where i, j, and k are unit vectors along x, y and z directions respectively. 


2.1.3 Energy Equation 

The first law of thermodynamics applied to a perfect gas passing 
infinitesimal, fixed control volume yields the following energy equation 

■^ + V-E,V = ~-V-q + />f •V + V-(lI,^ -v) 
dt dt 


where Et is the total energy per unit volume given by 


E, 


e -) H 

2 


(2.7) 


( 2 . 8 ) 


through an 


(2.9) 


( 2 . 10 ) 
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=u^ + v~ + 


( 2 . 11 ) 


and U , ^ is the stress tensor given by 


II„ =-M,; 


dll. du , 


+ - 


dx , dx. 


‘J 




ditk 

dx. 


(ij,k=l,2,3) (2.12) 


The first term on the left-hand side of the energy equation represents the rate of 
increase of total energy per unit volume in the control volume while the second term 
represents the rate of total energy lost by convection (per unit volume) through the 
control surface. 

The first term on the right-hand side of the equation is the rate of heat 
produced per unit volume by external agencies while the second term is the rate of 
heat lost by conduction (per unit volume) through the control surface. The third term 
represents the work done on the control volume (per unit volume) by the body forces 
while the fourth term represents the work done on the control volume (per unit 
volume) by the surface forces. 

For a Cartesian coordinate system, equation becomes 


^ + /z w) + ^ + pu - - vr^ - + q^) 

dt dt ox 

+ -^(e,v + pv- - vtyy - ) 

+ -|- (i:, w + pw- - wr,, + 0 


which is in conservation-law form. 

q = 9.> + ?yj + 9zk 
In terms of heat conduction 


(2.13) 


(2.14) 
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=-k 






q, =-k 


ydz , 


(2.15) 


2.1.4 Vector Form of Equations 

The governing equations for application of finite difference algorithm were taken to 
be two-dimensional unsteady, compressible Navier-Stokes equations without body 
forces, external heat addition or turbulence model Written in conservation form in 
Cartesian rectangular coordinates, these equations can be expressed as 


au aE ^ aF_Q 

dt dx dy 

where 


(2.16) 


P 

p u 
pv 

- 

pu 

pu^ +p-r^^ 
puv-v^y 

{E, +p)u-uT^^ -vTyy+q, 

pv 

pUV-ty^ 
pV^ +p-tyy 

(E, +p)V-UT^y -Vryy+qy 


(2.17) 


where 


E,=p 


e + - 


V 


2^ 
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2 du 

dx dy j 




2 9v 5 m 
dy dx 


= /“ 


du dv 
\^dy dx 


= T, 


and Sutherland’s law for viscosity coefficient 


H = II 


7 ; +110 

V^o y 


(2.18) 


r +110 

and the perfect gas relationship for calculation of pressure have been employed. The 
thermal conductivity was calculated from 


k =■ 


Pr 


(2.19) 


2.2 Boundary Conditions 

The boundary conditions specified for the computational domain are as follows 
1 . At the leading edge of the plate 


M = 0 

V = 0 

r=r 


( 2 . 20 ) 


2. At the inflow boundary (except the leading edge) and upper boundaries 


v = 0 

P = P. 
T = T 


( 2 . 21 ) 
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Figure 2.1: Application of Boundary Condition 


3. All the properties on the outflow boundary (except i = 0 and i = IMAX) 
are calculated based on an extrapolation from the two interior points, at 
the same i location. For example, u is determined as follows 

“(UMAX) = (2 “ (UMAX-i))~ “(i.JMAX-2) (2.22) 

4. On the lower wall, the no-slip condition (u = v = 0) is specified on velocity 
along with 


dn 


^0 


(2.23) 


and adiabatic wall condition 
8T 


dn 


0 


(2.24) 
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5. At the nozzle exit 


u =0 


P = Pje, 
T = T 

jet 


(2.25) 


2.3 Active Control Strategy for Transverse Injection 

The boundary conditions used for designing a controller employ the strategy of active 
control. An active control system will give the desired input as compared to a ‘passive 
control system’, where the output of the system may not be the desired one. The 
feedback control system used is shown in Fig. (2.2) 


M«> 1 



Figure 2.2: Schematic of Slot Injection Geometry with Control System 


1 





where Ae= — - (2.26) 

\P<o V P’a J 

is the error in the pressure ratio at a specific point and K is the controller gain. 

The control input to the plant is the change in total pressure that is calculated 

from 

-P; =AP,=K{Ae) (2.27) 


The new Pq and To (held constant) are used to calculate the properties of secondary jet 
stream. Assuming the flow in the nozzle to be isentropic 



where Pe = psep/2 [Ref. 4] (2.29) 

where pe is the exit pressure of nozzle and psep is the separation pressure, which is 
evaluated at a grid point just before the slot. 



(2.30) 

(2.31) 



(2.32) 


Tjet, Vjet, Pjet are the jet values at the nozzle exit. These are specified as the boundary 
conditions at the slot exit along with u = O.At the start of iteration pe and pjet will not 
be equal. They will be changing from iteration to iteration and becomes equal when 
the solution reaches the steady state. 
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Chapter 3 

Numerical Solution Procedure for Transverse Injection 

3.1 Numerical Solution of Navier-stokes Equations 

The finite difference method has been used to solve the governing flow equations 
presented in Chapter 2. The partial differential equations are converted to difference 
equations by discretization of governing equations on a finite grid. The flow field 
variables are evaluated at each node and after each time step. 

3.1.1 Discretization of Derivatives 

Each term of governing partial differential equations can be written in a finite 
difference form using taylor series expansion. For example, referring to Fig. (3.1) 
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Figure 3.1: Finite Difference Grid Representation 


^du'' 


\dyj 


(u) 






yM - y, 




’■i-uj 


yi - yi-i 


{Forward difference) 
{Backward difference) 


- {Central difference) 

yM - y.-i 


(3.1) 


3.1.2 Grid Generation 

Uniform and non-uniform grids are used to obtain the results. Fig. (3.2) shows the 
grid geometries used to get the results. Non-uniform grid in Fig. (3.2B) can be used to 
resolve the boundary layer more clearly by including a finer grid near the wall and a 




coarser grid near the upper boundary. The following transformation^*^ is used to get 
the stretching in the y direction. 




(a + l)-(a-l) 




va-ly 




(3.2) 


-ly 


+ 1 


Here ‘a’ is the parameter to check the fineness at the wall. As ‘a’ approaches ‘ 1 ’ more 
grid points are clustered near the surface. 

The grid in Fig. (3.2C) is used to resolve the near wall effects around the 
secondary jet. The following transformation^^ is used to get the stretching. 


Jc[;] = D 


^ ^ ((>(xl[/]-.4)) 

Sink ipA) 


(3.3) 


where 


.4 = — In 
lb 


l + (e* -l)f— 1 
\ ^ J 

1 + («■*-!( 


D'' 

L) 


(3.4) 


‘D’ is the length from the leading edge of the plate where clustering is required and 
b = 0 to 00 , is the parameter to check the fineness. When b = 0, no clustering will be 
enforced. L and H are the horizontal and vertical lengths of the computational domain 
respectively. 
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( C ) NON-UNIFORM GRID (12x12) 

Figure 3.2: Uniform and Non-Uniform Grid Geometry 


3-13 MacCormack Scheme 

The governing equation (2.16) is rewritten in vector notation as 

au aE aF .. .x 

dt dx dy 

Steady state solution to the above equation was obtained using the time-dependent 
explicit two-step predictor-corrector finite difference algorithm of MacCormack, 
which has been proved to be a dependable technique for producing numerical 

solutions to wide variety of complex fluid flow problems. 

The flow field variables are advanced at each grid point (i, j) in steps of time 

as 
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u; 


:U' + 

i,j 




ydt j 


At 


(3.6) 


where 


fdv) 

1! 


fau^ 

t 


r+At 

I at) 

= — 



+ 


2 

av 


[at J 


[at; 

'.7 


(3.7) 


is the average of derivative values of U of predictor and corrector steps. 

After each predictor or corrector step, the physical variables are obtained from 
vector U as follows 


P..J = 


‘.y 


‘.j 


",y 


'.y 


Ik 

U. 

Ik 

U. 

U4 

u 

u 


2 , 2 

4 “^y +^..y 


(3.8) 


Having obtained p, j, u, j, and the remaining flow field variables can be 
obtained as follows 


r = 

^ ij 

^/.y 




Cy 


Pu = 

Pu 

RT^ 



(t 

^ id 

Pij = 

0 

L 


Tq +110 

Tuj +no 


k 


i.j 


Pr 


(3.9) 
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3.1.4 Discretization of Boundary Conditions 

1 he X derivatives on the inflow boundary are obtained by forward difference and on 
the outflow boundary by backward difference. Similarly the y derivatives on the plate 
surface arc obtained by forward difference and that on the upper boundary by 
backward difference. For example, on plate surface 


^y. 


( 0 ./) 




(3.10) 




ym-y[0] 

The value of psep is calculated by extrapolation from the two interior points 

Ps.n= - P[2]U] (3.11) 


3.1.5 C Calculation of Time Step 

Because of the explicit nature of the procedure, a time step-size limitation exists for 
stability, 'fhe time step-size is calculated from Courant-Friedrichs-Lewy (CFL) 
criterion’^* and is given by 






At Ay Ay 


1 1 o • 

H :r 4- 2 V 


1 1 
+ - 


Ax^ Ay 


(3.12) 


where 


v' = max 


Pu 


(3.13) 


and the time step 


At = min[iiC^{(At CFL ) tj )j 


(3.14) 


for 0.5 <-K, ^ 0.8 .Kracts as a ‘fudge factor’, to make the solution stable. 
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3.2 Solution Procedure 


The steady state solution of the flat plate is taken as the initial condition. Boundary 
conditions for the jet are obtained from control strategy and MacCormack scheme is 
applied to calculate the flow field variables at the next time step. The procedure can 
be summarized as follows 

1 . Initialize the variables u, v, p and T at all the grid points obtained from the 
steady state solution. 

2. Find the other variables p, Et, e, p and k using the present value of u, v, p and T 

3. Calculate the time step by the CFL criterion. 

4. Calculatc the flow field variables after the predictor step using MacCormack 
scheme. 

5. Get the jot boundary conditions from the control strategy. 

6. Apply the boundary conditions. 

7. Repeat steps 4, 5 and 6 for the corrector step. 

8. Repeat the above steps until steady state is reached. 
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Chapter 4 


Results and Discussions 

I'he system used for computations is SG3400 workstation of IIT kanpur computer 
center. The program converges in 20 minutes for a grid of 70 x 70 for 8000 iterations. 

4.1 Validation of the code 

The qualitative validation of the code developed is done by comparing the computed 
surface pressure distribution with that shown in Ref. [13]. The computed surface 
pressure distribution in Fig. (4.1) shows a pressure rise because of separation shock 
followed by another pressure rise because of bow shock. Also the pressure fall 
because of separation downstream of jet and a pressure rise because of recompression 
shock can also be seen. The geometry selected for comparison consists of a slot of 
width W“0.0I89 L. The freestream values are taken to be T=288.16 K, p=101325 
N/ml The Reynolds number based on L is approximately 65500 and an uniform grid 
of 376x 1 82 has been taken to make the comparison 



4.2 Results of closed loop control study 

The configuration selected for control system study consists of a flat plate of length 
2L . The width of the nozzle at the slot exit is w=0.15L.The freestream values are 
taken as f— 288.16 K, p=101325 N/m^ This has been selected because it gives laminar 
flow over the entire flat plate with Reynolds number about 1000 so that computational 
running times are short as compared to high Reynolds number applications. The study 
has been performed for Mach numbers 3.5 to 4.5 The primary and secondary fluids 
are assumed to be air. 

The result in Figs. (4.2) to (4.5) shows the numerical experimentation done to 
design the controller for different desired pressure ratios. In Fig. (4.2), for a desired 
pressure ratio of 3.0, the controller values of 2.0 and 10.0 were not able to give the 
desired pressure ratio. A minimum value of K=30 is required to get the desired 
pressure ratio of 3.0. Similarly in Fig. (4.3), the controller values of 10.0 and 30.0 
were not able to give the desired pressure ratio of 4.0. A minimum value of K=50 is 
required to get the desired pressure ratio of 4.0. Likewise, a minimum value of K=100 
and K=200 are required for desired pressure ratios of 5.0 and 7.0 respectively and are 
shown in Mgs. (4.4) and (4.5). A value of K=100 can be the best one, to use for 
desired pressure ratios upto 5.0 for M=4.0, as more stagnation pressure will not be 
required in the transient state. A value of K=200 will be the best one for desired 
pressure ratios upto 7.0. But, the injection total pressure can be very high in the 
transient state for low desired pressure ratios. 

Figs. (4.6) to (4.9) shows the plots of desired output, the control input, 
injection Mach number and injection pressure variation with time for different values 
of K and M=4.0. As can be seen from the plots that with increase in the value of K, 
more injection stagnation pressure is required in the transient state, thus making the 
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flow to get choked for sometime, though the nozzle doesn’t need to be choked in the 
steady state. This is the advantage of using an active feedback control system, as for 
sonic cases of desired pressure ratios, it only requires an unchoked flow to get the 
desired picssure ratio. This is not the case with an open-loop control, used in all the 
references where they always used a choked flow. 

Fig. (4.10) shows the results for M=4.0 and K=200 for different desired 
pressure ratios. From Fig. (4.10b), for same ffeestream Mach number, in order to get a 
higher desired pressure ratio, the injection total pressure has to be increased. A higher 
injection total pressure means a greater secondary mass flow rate, thus increasing the 
strength ol bow shock wave, which gives a greater pressure rise downstream of the 
shock. As can be seen from Fig. (4.10c), the nozzle was not choked when desired 
pressure ratio was 4.0. But when the desired pressure ratio was 7.0 the nozzle gets 
choked, thus indicating a greater secondary mass flow required to get a higher desired 
pressure ratio. 'ITiough the nozzle was choked, increasing the injection total pressure 
can increase the mass flow rate. 

Fig. (4.11) shows the variation of different quantities with changes in Mach 
nurnher. ilicse plots can also be used to study the robustaess of this control system to 
changes in input quantities, in particular freestream Mach number, which is always 
possible in a physical situation. From the plots it can be said that the control system 
will give the desired steady state even when there are slight changes in input 
quantities. Though the control system was designed for M=4.0, it can sustain changes 
in Mach number from 3.5 to 4.5 for a desired pressure ratio of 4.0 and K=100. From 
Fig. (4. lie), wc can see that for the same desired pressure ratio, for a smaller 
freestream Mach number, we need a greater injection Mach number and stagnation 
pressure indicating a greater secondary mass flow rate. 


33 



Figs. (4.12) to (4.15) shows the Mach number contours for a transient and steady 
state for M— 4.0, I 4.0 and K=100. The separation regions before and after 
the slot can be seen from contours. Figs (4.14) and (4.15) shows the Mach number 
contours near the slot for a transient time and steady state. Fig. (4.15) shows the 
increase in the separation region as the solution reaches the steady state. Shear stress 
contours are shown in Fig. (4.16) and temperature contours in Fig. (4.17) respectively 
for M=4.0, Pj ! p„=^S) and K=100. Fig. (4.18) shows the surface pressure 
distribution for a uniform and non-uniform grid of 70x70. The non-uniformity is 
employed only in y-direction (a=2) with a uniform grid in x-direction. Increasing the 
fineness near the wall by employing a non-uniform grid shows that, the separation 
regions can be best captured. Also the non-uniform grid shows a higher pressure rise 
because of leading edge shock wave as more grid points are employed near the wall. 
Fig. (4.19) shows the comparison between a uniform and non-uniform grid of 70x70 
with non-uniformity employed only in x-direction. The plots shows that by employing 
more grid points before and after the slot, the separation regions can be captured more 
clearly. 

4.3 Identification of System Transfer Function 

The governing flow equations (Navier-Stokes) are of non-linear nature. We have to 
make some drastic approximations to make the equations linear about an equilibrium 
point in order to get the transfer function of the system. This is not a practical way of 
getting the transfer function of the system. Instead the modeling of such complex 
systems about an equilibrium point can be done by using techniques of modem 
control theory. This type of modeling can be done based on the input-output 
relationship of the system. The time histories of the output (desired pressure ratio) and 
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input (injection total pressure) are used to model the system by transforming them 
into the frequency domain and the resulting frequency responses are used to obtain 
the transfer function. Fast Fourier transform was used to transform the time history 
data of input and output to the frequency domain assuming the system to be linear. 

Assuming the differential equation associated with the system to be linear, the 
fast fourier transform of input and output was plotted for different Mach numbers and 
is shown in Fig. (4.20) for Pj / p„= 4.0 and K=100. Fast Fourier transform was 
obtained using MATLAB and the frequency vector used to get the fast fourier 
transform was taken from page 327 of Ashish Tewari The approximate transfer 
function of the system for M=4.0 was obtained as 


6’(.v) - 


3.1078 xlO'^ + 2.7x10^ + 3.2563 x lO'^ 
5^ +8.2901xl0’s + 1.8714xl0'® 


(4.1) 


that suits best in the frequency range 10* < o) < lO" .The comparison between the 
exact and approximate frequency response is shown in Fig. (4.21). This transfer 
function can be used to improve the system performance by designing a PID 


controller. 
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Figure 4.3: Variation of pressure ratio with time for pj 4.0 and M=4.0 
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(a) Pressure ratio 



(b) Injection total pressure 


Figure 4.6: Variation of different quantities with time for different K values and 
Pjl P«~ ’ M=4.0 
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(b) Injection total pressure 

Ingurc 4.7: Variation of different quantities with time for different K values and 
Pj Pm =4.0 ,M=4.0 
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Figure 4.12: Mach number contours at a transient time (250 iterations) for M=4.0 
p,/p^=4.0 andK=100 
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Figure 4.14: Mach number contours near the slot at a transient time (100 iterations) 
for M=4.0 Pj / p^= 4.0 and K=100 


53 



l’'igure 4. 1 5; Mach number contours near the slot at steady state (7000 iterations) for 
M--4.0, p,, / = 4.0 and K=100. 
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Figure 4.16: Shear stress x^y contours for M-4.0, p, / 


4.0,K=100 






Figure 4.1 8: Comparison of surface pressure distribution for a Uniform and 
Non-uniform grid (70x70) 


57 




Figure 4.19: Comparison of surface pressure distribution for a Uniform and 
Non-uniform grid (70x70) 
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Figure 4.21 : Comparison of Fast Fourier and Approximate Transfer function 
frequency response for M=4.0 Pi! ^.0 , K=100 










Chapter 5 


Conclusions and Suggestions for Future work 

5.1 Conclusions 

A numerical study to design a linear feedback control system for the case of 
transverse injection over a flat plate in supersonic flow that gives a desired static 
pressure ratio upstream of the slot has been carried out by integrating the two- 
dimensional unsteady Navier-Stokes equations. Isentropic flow was assumed in the 
nozzle. The computed result was compared with a reference to check for its 
validation. The study has been performed for different desired static pressure ratios 
and controller constants. Results depicted that for same freestream Mach number, a 
higher desired pressure ratio requires greater secondary mass flow, which increases 
the strength of bow shock and thus giving a higher desired static pressure ratio. From 
the results, we can see that the nozzle doesn’t need to be choked for some cases of 
desired pressure ratios thus showing the advantage of using an active feedback control 
system, as it requires an unchoked flow to get the desired pressure ratio. This is not 
the case with an open-loop control, used in all references where they always used a 
choked flow. Also such a control system gives the exact desired value even for small 
changes in freestream Mach numbers thus showing the robustness of system to small 
changes in freestream Mach numbers. The transfer function of the system has been 
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obtained for a particular case which can be used to improve the performance of 
control system by including a non-linear controller. 

5.2 Future Scope of the Work 

Much of scope exists for extending this work. The control system designed finds its 
application in scramjet engine where there is considerable importance of the 
maximum pressure for mixing enhancements of fuel and air and also in thmst vector 
regulation. The same concept of controlling pressure can be extended to control heat 
transfers for cooling the surfaces and controlling skin friction for drag reduction in re- 
entry vehicles. 

The controller design in this thesis was carried out by numerical 
experimentation. But, the controller design can also be carried out by transfer function 
approach. Once the transfer function of the system is obtained about an equilibrium 
point, the same can be used to design the controller, which can be an extension of this 
thesis. The code can also be extended to include the turbulence modeling for the case 
of turbulent flows and also can be extended for a three-dimensional problem, by 
including the third spatial parameter. The code can be made more accurate by 
including non-equilibrium effects when solving the problem in a higher Mach number 
range. 
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